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The present work deals with a set of problems in isotope shifts of neutral barium spectral lines. 
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mean-square radii with our electronic factors for the 6s5(i ^Di ^2 — 6s6p ^P° transitions produces 
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I. INTRODUCTION 


Q. 


who 


The first isotope shift (IS) measurements on barium have been done by Arroe 
studied the resonance transition 6s^ ^So — 6s6p of neutral barium located at a wavelength 
of Ai = 553.7 nm. One of the first attempts to obtain the differences of nuclear mean-square 
charge radii of the radioactive barium nuclides was done by Fischer et ai who 

investigated experimentally and theoretically the isotope shift in the Ba II resonance line 
6s ^S'i/2 — Qp ‘^Pi/2 at A = 493.4 nm. Their results are often used by other authors. Later, 
the Doppler-free spectroscopy was explored by Nowicki et al. js, 4| who studied the neutron- 


dehcient isotopes ^^^Ba and ^^^Ba. Using the same procedure Bekk et al. pursued the work 
of Nowicki et al. with other unstable isotopes. Concurrently but independently, Baird et al. 
B proposed a s.m.lar exper.meat. Years later, thanks to h.gh-reso.ntion laser speetroseopy 
Grundevik et al. [7[ explored the far-red transitions between the 6s5d and 6p5d configurations 
and, in continuity, investigated the spin-forbidden transition 6s^ ^-S'o — 6s6p ^P° at A 2 = 
791.3 nm js]. In a work of Mueller et al. jo], collinear laser spectroscopy has been connected 
to the mass separator ISOLDE-II at CERN. One of the aim of this experiment was to 
extend the knowledge on fundamental nuclear properties into regions far from stability. In 
that context, isotope shifts of barium isotopes in the mass range 122-146 have been measured 
for the Al atomic transition. Finally, Wijngaarden and Li 10| (re)measured IS of the same 


transition using a ring dye laser and obtained the most recent and precise IS values o 
line. Many other measurements were also reported on highly excited states of Ba I [iT 


this 

12 |. 


The theoretical barium studies are far less advanced. In 1974, Trefftz et al. 13| performed 


calculations on various states of barium using MCHF wave functions generated with the 


code of Froese Fischer [l^ but did not study IS. King and Wilson [15|] used a modihed 
version of the latter program that includes the mass velocity and Darwin terms in order to 
calculate electron densities at the nucleus. Besides, Fricke et al. 161 and Olsson et al. [l^ 
used MCDHF wave functions to compute electronic F factors. Finally, in an unpublished 


work, Kozlov and Korol 


18| calculated the field and mass shifts of Ba I and Ba II using 


second order many-body perturbation theory to take core-valence and valence correlations 
into account. 


Our interest has been picked up by the paper of Dammalapati et al. [l9| who reported 
the first laser spectroscopy measurements of the 6s5d ^Di 2 — Qs6p ^P° transitions for several 
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isotopes. Observing that the resulting isotope shifts strongly deviate from their expected 
behavior for odd isotopes in an analysis based on King-plots, the authors pointed out that 
there were no theoretical calculations available for comparison. 


II. ISOTOPE SHIFT THEORY 


The main ideas of the isotope shift theory are hereafter outlined. The interested reader 
should look at the pioneer works of Shabaev 2Cll-l24| and Palmer 25(] who expressed the 


theory of the relativistic mass shift as used in the present work, 
the relativistic recoil operator was derived by Gaidamauskas et al. 


he tensorial form of 


26| . Based on those 


developments a module, called Ris3 (Relativistic Isotope Shift) was designed |27| for the 
revised version ofGRASP2K package 


A. Mass shift 


The hnite mass of the nucleus gives rise to a recoil effect, called the mass shift (MS). The 
nuclear recoil corrections within the (aZ)^m^/M approximation are obtained by evaluating 
the expectation values of the operator 

1 A / aZ ( (cKi-rdrA \ 

L 1 Pr P. - — 


Hms — 


2M 


Pi ■ Pi 


OLi 


■Pi 




( 1 ) 


Separating the one-body (i = j) and two-body {i ^ j) terms that respectively constitute the 
normal mass shift (NMS) and specihc mass shift (SMS) contributions, the Hamiltonian (fTj) 
becomes 


Hus — h^NHS + Hsus ■ 

(2) 

The (mass-independent) normal mass shift iPxMS and specihc 

are dehned by the following expressions 

mass shift iPsMS parameters 

K]SSMS = dT(T ipNMS d') , 

(3) 

-^SMS = dT(T i/sMS d') • 

(4) 


When discussing a transition isotope shift, one needs to consider the variation of the mass 
parameter from one level to another. The corresponding line frequency isotope mass shift is 
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written as the sum of the NMS and SMS contributions: 


r. A A' r. A A' ^ A A' 


( 5 ) 


with 


(M'-M\AKms 
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MM' ) 


AKms. , 


( 6 ) 


where AJ-^ms = (A'u^ms - A^^^ms) is the difference of the Kms (= A'nms + A'sms) param¬ 
eters of the levels involved in the transition k. For K, the unit (GHz u) is often used 
in the literature. As far as the conversion factors are concerned, we use AKMs[^eEh\ = 
3609.4824 A(Ams[GHz u]. 


B. Field shift 


The energy shift arising from the difference in nuclear charge distributions between two 
isotopes A and A' for levels i = {i, u) involved in transition k, the frequency held shift (FS) 


of the spectral line k can be written as 


29 


3l| 


A,A' ^-^n,FS “ 

=- T, 


e,FS 


Ek 5{rY’^' • 


(7) 


Ek is the line electronic factor 

proportional to the change of the total probability density at the origin associated with the 
electronic transition between levels i and u. In this approximation, the hrst-order frequency 
held shift becomes 


r A,A' XI 2\A,A' 

^^fc,FS — Ek S{r ) 


z 
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(9) 


C. The total isotope shift 


It is easy to estimate the total line frequency shift by merely adding the mass shift ([S]) 
and held shift ([7]) or ([9]) contributions 


6u 


A,A' 
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'M' - M" 
MM' . 


AKms + Fk 6{r^) 


2\A,A' 


III. COMPUTATIONAL PROCEDURE 


The multiconfiguration Dirac-Hartree-Fock (MCDHF) method 3^, as implemented in 


the program package grasp2K 


28 


33| . is employed to obtain approximate wave functions 


that we will refer to as atomic state wave functions (ASFs). An ASF is represented by a 
linear combination of configuration state functions (CSFs) with same parity P, total angular 
momentum J and component Mj along 2 ;-direction 


NCSFs 

T(7PJMj)= ^ c,-$(7 ,PJMj), 
i=i 


( 12 ) 


where {cj} are the mixing coefficients and { 7 . 4 } the sets of quantum numbers needed for 
specifying unambiguously CSFs. The latter are built from single-electron orbital wave func¬ 
tions. Applying the variational principle, the mixing coefficients and single-electron orbital 
wave functions are optimized simultaneously in the relativistic self-consistent field (RSCF) 
method. The energy functional is estimated from the expectation value of the Dirac-Coulomb 
(DC) Hamiltonian 341 


N N 

Hoc = X) (c ■ Pi + (A - 1)A + V (r*)) -F 


2 = 1 


i<j Dj 


(13) 


where V (rj) is the monopole part of the electron-nucleus interaction, a and f3 are the (4 x 4) 
Dirac matrices and c is the speed of light (c = 1/a in atomic units, where a is the fine- 
structure constant). 

To effectively capture major electron correlation, CSFs of a particular parity and sym¬ 
metry are generated through substitutions within an active set of orbitals occupied in the 
reference configuration. As regards of the hardware and software limitations, it is obviously 
impossible to include all CSF in the sense of a CAS expansion. The CSFs expansions have 
to be constrained so that the major correlation excitations are taken into account. In our 
calculations an approach based on single (S) and restricted double (rD) substitutions was 


applied 


. Restricted double substitutions limit the excitations to maximum one hole in 


the core. In the case of barium ([Xe] 6 s^, [Xe] 6s6p, [Xe] 6 p^, [Xe] 6s5d), double excita¬ 
tions are applied to the valence shells but with the restriction that at most one electron is 
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TABLE I. Sample from several sources of the nuclear rms radii and difference of the nuclear mean- 


square charge radii relative to 

in fm and fm^, respectively. 




Semi-empirical formula (11411 

From Angeli [3^ 

From Angeli and Marinova [3^ 

(,2)1/2 

j(,2)138A 

(,2)1/2 

^(,2)138A 

(,2)1/2 

^(,2)138.A 


4.8901 

0.0 

4.8385(48) 

0.0 

4.8378(46) 

0.0 


4.8798 

0.1006 

4.8326(48) 

0.0571 

4.8314(47) 

0.0609(2) 

i36Ba 

4.8692 

0.2040 

4.8327(48) 

0.0561 

4.8334(46) 

0.0422(2) 

i35Ba 

4.8586 

0.3071 

4.8273(48) 

0.1082 

4.8294(47) 

0.0812(3) 



S (r 

2 ) 138,A from different sources 



Ref. \£ 

Ref. 

Ref. 

Ref. [37] 

Ref. [7] 

Ref. [10] 

i38Ba 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

i37Ba 

0.049(2) 

0.059(4) 

0.067 

0.072(6) 

0.059(6) 

0.049 

i36Ba 

0.034(3) 

0.041(8) 

0.061 

0.044(4) 

0.042(7) 

0.034 

i35Ba 

0.066(5) 

0.079(12) 

0.108 

0.088(7) 

0.080(10) 

0.065 

i34Ba 

0.044(7) 

0.053(16) 

0.095 

0.051(7) 

0.056(10) 

/ 


substituted from Is2sp3spd4:spd5sp [Xe] shells, the other(s) involving the 6s, 5d or 6p va¬ 
lence shells of the considered reference conhguration. All occupied orbitals in the reference 
conhguration are treated as spectroscopic and are obtained in the DHF single-conhguration 
approximation. These orbitals are frozen in all subsequent calculations. The J-levels be¬ 
longing to a given term were optimized simultaneously with standard weights through the 
EOT scheme and the set of virtual orbitals is increased layer by layer. Starting from the 
n = 9 correlation layer, the core is fully opened but the new CSF are generated by single 
excitations from the reference. The effect of adding the Breit interaction to the DC Hamil¬ 
tonian flTdl) is estimated to be much smaller than the uncertainty in the transition isotope 
shift parameters with respect to the correlation model. This interaction has been therefore 
neglected. 
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TABLE II. Experimental measurements of isotope shifts for the transitions — 6s6p and 


6s 2 1^0 - 6s6p (in GHz). 



138 

137 

136 

135 

134 

6s^ — 6s6p ^Pi 

Ref. [1^ (1995) 

0.0 

-0.21515(16) 

-0.12802(39) 

-0.25929(17) 

/ 

Ref. (1979) 

0.0 

-0.2150(7) 

-0.1289(5) 

-0.2609(7) 

-0.1430(5) 

Ref. (1979) 

0.0 

-0.2147(5) 

-0.1275(13) 

-0.2587(7) 

-0.1428(12) 

Ref. [1] (1950) 

0.0 

-0.16(2) 

-0.066(20) 

-0.22(2) 

-0.13(2) 

6s^ ^S'o — 6s6p ^Pi 

Ref. [8] (1983) 

0.0 

-0.1834(10) 

-0.1092(10) 

-0.2199(10) 

-0.1223(25) 


IV. ISOTOPE SHIFTS OF NEUTRAL BARIUM 


A. A large “choice” of nuclear radii 


The reliability of the FS values obtained with the ab initio electronic F factor flHl) is a 
function of the accuracy of the calculations but also of the level of confidence on the nuclear 
data Table H reports, nuclear root-mean-square (rms) charge radii and mean- 

square charge radii differences from several sources, taking ^^®Ba as the reference isotope. 
The aim of this table is to illustrate the difficulties to reduce uncertainties on the nuclear 
rms charge radii. The first set of columns is obtained with the semi-empirical formula 

Rrms = = 0.836A^/3 + 0.570 fm if A > 9, (14) 


and is compared with the values compiled by Angeli 35j and Angeli and Marinova 36j . The 
positive sign of 6 (r^) (according to the convention = (r^)^ — with A > A') 

indicates that the neutron deficient isotopes of barium have smaller rms radii than ^^®Ba. 
The second-half of the table reports the values deduced from different tech¬ 


niques. Bekk et al. measured IS by laser induced resonance fluorescence on an atomic 


beam of a sample of barium. In continuity of the work of Nowicki et al. { 4 ], they connect^ 


^P° transition in Ba I to the data of Fischer et al 


■ s 


their IS measurement of the ^Sq 
on the A = 493.4 nm transition of Ba II via a King plot. Thanks to the electronic factor F 
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and SMS value deduced in 


j^, Bekk et al. 


were able to get their own F factor. They 


deduced the NMS through the relation 




nie 

M' 


me 

M 


'exp • 


(15) 


and, assuming that 5 i^sms = (0±l)(5r'NMS j^l, obtained the nuclear mean-square charge radii 
values. As pointed out by the authors themselve^ the latter approximation is the major 
source of the uncertainties on 6{r‘^). Baird et al. ^ combined their own results with other 
optical data and electronic x-ray measurements. Muonic x-ray measurements of the nuclear 
charge radii and 6{r‘^) values were presented by Shera et al. [37|. By comparing IS data on 
the ^S'o — transition with the muonic 6{r^) values, Shera et al. extracted the electronic 
factor. Grundevik et al. j?] reevaluated the SMS values in the Ba I and Ba II resonance 


and Li 


transitions and obtained results strongly consistent with On their side, van Wijngaarden 
10| used the procedure proposed by jd] on their new IS measurements. Using the 


electronic factor of Fischer et al. |2|, they deduce a new value of the variation of the nuclear 

n 

mean-square 6 (r^), in very good agreement with 

It is worthwhile to notice that in some cases the rms charge radius decreases when 

the number of neutrons increases along the isotope chain. This will never be reflected when 
using the semi-empirical formula flTdl) whose values increase monotonically. Therefore, in 
present paper, this approach will be left aside when rms charge radii are needed, in favor 


of data coming from the literature such as Angeli (and Marinova)’s compilations 35 


In addition, table [I] demonstrates the difficulties to isolate the nuclear rms radius and sheds 
light on the remaining large uncertainties on the 6{r‘^) nuclear data for these systems. 


B. Some well-known transitions 

The resonance transition 6s^ ^S'q —6 s6p ^P° is maybe the most well-known in barium and 
is, together with the intercombination line 6s^ — Qs6p a good starting point of our 

analysis. ^^®Ba is the most abundant isotope of barium on Earth with its 82 neutrons and as 
such is often chosen as the pivot; the IS relative to the isotope 138 are given from A = 137 
to 134 in table m Some details behind those measurements are given in the introduction of 
the present paper. The values of Arroe [l|, originally given in cm“^, were converted so that 
the given error bars of ±0.7 x 10“^cm“^ become ±21 MHz. 













TABLE III. Experimental measurements of neutron-rich isotope shifts for the transitions — 

6s6p (in GHz) and their related difference nuclear mean-square radii with the reference isotope 


i38Ba. 



141 

140 

139 

6s2 1^0 - 6s6p 




Ref. [39] 

1.505(8) 

1.075(6) 

0.473(6) 

Ref. [9] 

1.505(5) 

1.075(3) 

0.473(3) 


^^^2^138,141 

^^^2^138,140 

^^^2^138,139 

Ref. [3^ (1981) 

-0.395(13) 

-0.281(9) 

-0.124(5) 

Ref. (1983) 

-0.395 

-0.281 

-0.124 

Ref. [4^ (1992) 

-0.440(1)(13)(1) 

-0.314(1)(9)(1) 

-0.1381(5)(41)(9) 

Ref. [36] (2013) 

-0.410 

-0.292 

-0.129 


The negative signs indicate that the isotope ^^®Ba has the lowest line frequency for each 
i38A'ga isotope pair (A' < 138). Assuming the dominance of the FS, these isotopes behave 
in the most current way, considering the density reduction with electronic excitations (i.e. 
(fz/pg"^ < 0 with A > A'). The consistency between experiments is rather good. 


In Neugart et al. 


and Mueller et al 


■ B. 


IS measurements are reported for neuron-rich 


isotopes 389Ba to ^^^Ba and the mean-square nuclear charge radii 6{r^) have been estimated, 


following the procedure of Bekk et al. in which the unknown SMS is taken to be of 
the order of the normal mass shift. These results are presented in the separated table HTTl 


The rela 


nerman 


dvistic coupled-cluster approach has been used by Martensson-Pendrill and Yn- 


4 of to calculate an electronic F factor that allowed the authors to revise values 
of nuclear charge radii using optical isotope shifts for muonic results for the 

stable isotopes Muonic results S{r‘^)^ were used to derive a Ksms parameter, itself 

used to extract (5(r^)opt. The error bars reflect the uncertainty in the SMS, in F and in the 
optical data. The agreement within the three sets of of 6{r‘^) is quite satisfactory. 
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1. Mass shift calculation 


The results of our calculations for the mass factors of transitions ^Sq — 6s6p 
and 6s^ — Qs6p are reported in table HVl For each transition, the hrst column gives 

the value of the i^MS parameters in GHz u, the second column is the value of the MS in 
MHz for the isotopic pair ^ multiplying the Kms parameter by the mass factor 

(1/Mi 38 — l/Mise)). Nuclear masses (M) are calculated by substracting the mass of the 
electrons and the binding energy from the atomic mass (Matom), using the formula: 


M(Ai Z) — Matom(^, Z) — ZrUe + Bel{Z) 


( 16 ) 


where the total electronic binding energy (expressed in eV) is estimated using 1^| 


^—6 


Bei{Z) = 14.4381 + 1.55468 • 10"*’ Z 


(17) 


Atomic masses are provided by 


43| 


TABLE IV. Values of the AAms parameters for the — 6s6p ^Pf and 6s^ ^5o — 6s6p ^Pf 

transitions (in GHz u) and values of the MS (in MHz) for the I38d36gg^ pair. 



6s^ 

AAms 

— 6s6p 

Ipo 

MS 

6s^ 

AATms 

— 6sQp 

3po 

MS 

DHF 

296.65 


-31.68 

242.69 

— 

-25.92 

6SrD 

-393.26 


42.00 

-125.10 

—)■ 

13.36 

7SrD 

-21.81 

—)■ 

2.33 

-203.07 

—)• 

21.69 

8SrD 

37.28 


-3.98 

-183.89 

—)• 

19.64 

9S 

61.66 


-6.58 

-184.19 

—)■ 

19.67 

lOS 

-59.65 

—>• 

6.37 

-180.10 

—)■ 

19.23 

ns 

-61.22 


6.54 

-179.14 

—)■ 

19.13 


Clearly, correlation effects play a major role and core excitations are non-negligible. How¬ 
ever, a comparison of the MS with the isotope shift values of the table [TTl reveals that the 
contribution of MS is really small and represents only 5% of the total isotope shift value for 
the resonance transition but reaches around 20% for the transition ^Sq — ^Pf. In conse¬ 
quence, for both lines one expects a much larger contribution of the FS. Furthermore, one 
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observes a much larger stability and better convergence of the transition parameter AJ’^ms 
for the spin-forbidden transition (with a 0.5% difference between the IQS' and IIS' values) 
than for the resonance transition suffering from oscillations, even in the largest calculations. 
Actually, the values of A'ms of the and ^5*0 levels are so close that the slightest change 
on the level parameter strongly affects the transition parameter APms- The accuracy of the 
latter is therefore hard to evaluate but its reliability is discussed further in relation with the 
5(r^) values that can be deduced from the experimental IS (see table |VT|). Behind the value 
of the mass parameter APms = —61.22 GHz u for the transition Ai, hides the sum of the 
APnms = 362.51 GHz u and AA'sms = —423.74 GHz u. By looking at them, it seems that 
the approximation Susus = (0 ± l)(5r'NMS proposed by Nowicki et al. j4| is not senseless. 


2. Field shift calculation 


The level field shift in Ba I is around 10“'^Ph for both states 6s^ ^S'o and 6s6p while 
the transition FS is 10“®Ph; a good accuracy is not easy to reach especially for a total binding 
energy around —SOOOPh- With that respect, the formalism (|9]) is more reliable in view of the 
extreme difficulty to obtain highly converged total energies. Furthermore, the perturbative 
approach offers the freedom to explore and seek for the best nuclear mean-square radius. 

Table |V] gives a chronological list of the experiments and calculations performed so far to 
determine the electronic P factor of the 6s^ ^S'q-6 s6p ^P° and 6s^ ^S'q-6 s6p ^P° transitions. 
In the second-half of the table are presented the results of our calculations, starting from 
Dirac-Hartree-Fock to our most correlated model. 


To our knowledge, it is the hrst time that such a highly correlated model is used on 
barium. This could be an explanation for the fact that our values present some differences 
with the literature. For the ^S'o — ^P° tranship our values are in favor of a large P factor, 
in line with Bekk et al. and Kunold et al. 441. For the spin-forbidden line, our P factor 
values are larger (13%) than any others available in the literature. The held-shift parameters 
presented in table |V] have converged within 0.1% for both transitions. 
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TABLE V. Comparison of different theoretical and experimental determinations of the electronic 
F factor. 



F (GHz/fnF) 


^Sq — 6s6p 

^So - 6s6p 

Theory 



Ref. [8j 

-2.34 

-2.34 

Ref. [1^ 

-2.99 

-2.55 

Ref. \M 

-2.996 

-2.546 

Ref. [1^ 



Ref. [IT] 

-2.998 

-2.544 

Ref. [40j 

-3.39(10) 


Experiment 



Ref. 

-3.929 


Ref. [37j 

-3.04(26) 


Ref. [1^ 


-2.59(22) 

Ref. [44j 

-3.99(65) 


This work 

GRASP2K 


(DHF) 

-3.48 

-2.76 

(lOS) 

-3.95 

-3.49 

(US) 

-3.95 

-3.49 


3. extraction 

Thanks to the AKus parameters of table HVl the electronic F factor of table |V] and 
the formula m, we propose to isolate the values of 6{r^) that would reproduce the total 
measured isotope shift value (see tables UTI and IIIID . For the transition 6s^ ^S'q —6 s6p ^P°, one 
uses the IS values of van Wijngaarden and Li that appear to be the most precise experimental 
values. It is possible to double-check the consistency of our results by extracting the 6{r‘^) 
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values from the experimental IS of the other transition, i.e. 6 s^ ^S'o — 6s6p ^P° [7|. These 
results are presented in table IVTI and look really promising. Error bars on reflect 

the uncertainty on the IS measurements. The ( 5 ^ 7 - 2 ^ 138,137 obtained from these two transitions 
differ by less than 1%. The (5^7-2^138,136 138,135 values resulting from the two lines 

agree within 8 to 5%, respectively. However, the ( 5 ^ 7 - 2 ^ 138,134 reveal a discrepancy of 

17%. 


TABLE VI. Extraction of 5{r‘^) values (in fm^) using experimental IS values (in MHz) from the 


literature and our electronic factors MS (in MHz) and F (in GHz fm^): (5(r^)=(IS—MS)/F. 



6s 2 1^0 - 6s6p ipf 

( ISexp 

MS ) 

IF 

-> 


i4iBa 

1505(8) 

-9.48 

-3.95 

-> 

-0.3832(20) 

i40Ba 

1075(6) 

-6.37 

-3.95 


-0.2736(15) 

i39Ba 

473(6) 

-3.21 

-3.95 

-> 

-0.1205(15) 

i37Ba 

-215.15(16) 

3.24 

-3.95 

->• 

0.0553(0) 

i36Ba 

-128.02(39) 

6.54 

-3.95 

-> 

0.0341(1) 

i35Ba 

-259.29(17) 

9.88 

-3.95 

-> 

0.0681(0) 

i34Ba 

-143.0(5) 

13.3 

-3.95 

-> 

0.0395(13) 


6s^ ^5o — 6s6p ^P° 






( IS 

MS ) 

IF 

-> 


i37Ba 

-183.4(1.0) 

9.49 

-3.49 

-> 

0.0553(3) 

i36Ba 

-109.2(1.0) 

19.13 

-3.49 

-> 

0.0368(3) 

i35Ba 

-219.9(1.0) 

28.90 

-3.49 

-> 

0.0714(3) 

i34Ba 

-122.3(2.5) 

38.83 

-3.49 

-> 

0.0462(7) 


6s 2 1^0 - 6p2 Spp 






( IS 

MS ) 

IF 

—> 

^^.^2^138,A 

i37Ba 

-331.7(5.0) 

21.49 

-7.09 

-> 

0.0498(7) 

i36Ba 

-199.0(3.0) 

43.33 

-7.09 

-> 

0.0342(4) 

i35Ba 

-396.1(5.9) 

65.45 

-7.09 

-> 

0.0651(8) 

i34Ba 

-219.0(9.9) 

88.94 

-7.09 

-> 

0.0433(14) 
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The isotope shifts of the 6 s^ ^5'n — ^Pq transition reported in table IVTl are taken from 


the work of Jitschin and Meisel 


45| . They resolved the IS for several highly excited states, 


using Doppler-free two-photon laser spectroscopy. They needed the relevant electronic F 
factors to extract the 6{r‘^) values and the only response was from Olsson et al. 0 . Based on 
the correlation model detailed in section HU we calculated the IS parameters of the ^Pq 
state. The errors bars are sensitively larger than for the two other experiments. However, 
the consistency between the three independent sets give us confidence in the reliability of 
our electronic parameters {Kus and F), which is the original point of this section. 

Let us refer to the table [I] in order to compare these newly extracted values with the 
available ones. The 6{r‘^) values from the semi-empirical formula as well as results of Baird et 
al. jo I are both out of range but there are some nice agreements with the other experiments. 
The results of Bekk et al. ^ and Grundevik et al. seem to be confirmed, especially 

for (5^7-2^138,137 ^ 0.055 fm^ but not for ,5^7-2^i38,i36^ which, our results are closer to 

van Wijngaarden and Li’s value ("^^^ 2 ^ 138,136 _ q Q 34 fni2). The , 5 ^ 7 - 2 ^ 138,135 jg 
disputable, but let us just underline that our result is closer to van Wijngaarden and Li as well 
as the value obtained by Martensson-Pendrill and Ynnerman ,5^7-2^138,135 _ o.0728(15)(22)(2) 
fm^ 40j. The discrepancy found for , 5 ^ 7 - 2 ^ 138,134 forbic^g us to draw any further conclusions. 
All our results for neutron-rich isotopes confirm the deductions of Neugart et al. 39|] and 
Mueller et al. jo] within 3% (see table m- 


V. THE ISSUES OF THE ^Pf - TRANSITIONS 


The level scheme of barium also exhibits low lying 6s5d ^^ 2 , ^Di and ^£>2 metastable 
states. However, the transitions 6s5d ^Di — 6s6p ^P° and 6s5d ^D 2 — 6s6p are far less 
studied and the first measurements of their IS have been performed by Dammalapati et 


al. 


19| . The latter are displayed in table rvTTTI using the ^^®Ba as the pivot. These authors 


also measured shifts of other pairs, reporting = —84.8(8) MHz and —80(1) MHz for 

the transitions — ^Di and ^P° — ^^ 2 , respectively, and = —75.3(5) MHz for 

the resonance line frequency. 

Following the computational procedure described in the section IIIIl the P and Pms 
parameters have been calculated and are reported in table I VI II 

It firstly appears that the electronic P factor is much smaller and the APms param- 
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TABLE VII. Values of the electronic F factor (in GHz/fm^) and of the AKms parameters (in 
GHz u) for the 6s5d ^Di — 6s6p and 6s5d ^D 2 — 6s6p transitions. 



F (GHz/fm^) 

AAms (GHz u) 

3Di - ipf 

3^2- ^Pf 

3_Di - Ipf 

b 

to 

1 

^0 

(DHF) 

0.31 

0.31 

937.12 

962.56 

(lOS) 

0.75 

0.74 

-284.34 

-324.09 

(US) 

0.76 

0.74 

-311.99 

-294.94 


eter larger than for the three transitions considered above (see table lYD- The mass shift 
value is = 33.32 MHz while the frequency FS is = 25.73 MHz (using 

5(r^)i38>i36_g_g3q fm^). Surprisingly, they are both of the same order of magnitude. The 
observation that FS strongly dominates the MS in high-Z systems cannot be obviously ap¬ 
plied to neutral barium. Moreover, we would recommend the greatest caution about the 
transition considered. Furthermore, the difference between the two most elaborated calcu¬ 
lations of the values presented in table IVIII is around 9% for the two ATFms parameters and 
0.5% for both electronic F factors. Thanks to those two observations, we could estimate 
errors bars of 5% on the IS values obtained from our electronic factors calculations for the 
6s5d ^Di — 6s6p ^P° and 6s5d ^D 2 — 6s6p ^P° transitions. 

Table IVHII compiles the results of Dammalapati et al. with our predictions using the 
(5(r^)i38A Qf other authors. The underlined data are obtained with the (5(r^) values that 
look the most reliable on the basis of the 3.ipo _ 15 '^ transitions. Table rvTTTI brings much 
interesting information. Most of our values are in complete contradiction with the results of 
Dammalapati et al. The pair ^38>i37gg^ jg consistent but only for the ^D 2 — transition. 
It is true that as regards to the proximity of the two levels 6s5d ^Di 2 , h is rather strange 
that the IS MHz for one transition is reduced by a factor 2 for the second 

one. All the and values show signs opposite to ours and that, whatever 

the (5(r^)^38A used. However, are consistent on the absolute scale. As for the 

its determination suffers from our incapacity to discriminate between right and 
wrong ^(7-2^138,134 values. 

In order to circumvent the problem, table IIXI presents the AKms values that would 
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TABLE VIII. Experimental measurements of isotope shifts for the transitions 6s5d ^Di —6s6p 
and 6s5d ^D 2 — 6s6p (in MHz) from Q compared with the IS calculated with the mean-square 
radii differences available in the literature and our electronic factors. 



138 

137 

136 

135 

134 

6s5d ^Di — 6s6p 






Dammalapati et al. [1^ 

0.0 

114(4) 

-59.3(6) 

39(4) 

-144.1(10) 

With the 6{r‘^) of: 






Baird et al. [6] 

0.0 

67.2 

79.5 

132.1 

139.5 

Bekk et al. 


61.2 

64.4 

110.1 

107.7 

Shera et al. 


71.0 

66.6 

116.9 

106.2 

Grundevik et al. 


61.2 

65.1 

110.9 

110.0 

van Wijngaarden et al. [1^ 


53.6 

59.1 

99.5 


6s5d ^D 2 — 6s6p 






Dammalapati et al. [1^ 

0.0 

69(3) 

-63(1) 


-143(1) 

With the 6{r‘^) of: 






Baird et al. [6] 

0.0 

65.0 

76.5 

127.2 

134.0 

Bekk et al. 


59.1 

61.7 

105.8 

103.0 

Sherra et al. 


68.7 

63.9 

112.5 

101.5 

Grundevik et al. [7] 


59.1 

62.5 

106.6 

105.2 

van Wijngaarden et al. [1^ 


51.8 

56.6 

95.5 



reproduce the experimental results of Dammalapati et al. [lOj]. In order to do so, we trusted 
the values of our F factors and adopted the most trustable 5{r‘^) and the nuclear masses from 
the literature. This is done with the help of the formula AiLMS= (IS—P(5(r^))/(1/Mi38 — 
1/M'). For the comparison, we deduced a value for ( 5 ( 7 - 2 ^ 136,134 making the difference 

(f(r^)136d34 _ ^(j,2^138,134 _ 138,136 _ Q QQ5 {yS?^ g^g fgj- ^(,,2^137,135 _ ^(^,2^138,135 _ 

(5(r2)i38,i37 _ q_qx 7 fm^. 
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TABLE IX. Extraction of the AA^ms that would be necessary to reproduce the experimental results 
of Dammalapati et al. [l^ . combining our theoretical F factor with the 5{r‘^) and nuclear masses 
from the literature. 


AA'ms = - F5{r^))/{l/M - 1/M'). 




(MHz) 

F 

(GHz/fm2) 

(fm^) 

(1/M - 1/M') 
(u-i) 

AAms 

(GHz u) 

6s5d ^Di 

— 6s6p 







138-137Ba 

114(4) 

0.76 

0.055 

-0.000053 ^ 

-1366.4 


138-136Ba 

- 59.3(6) 

0.76 

0.035 

-0.000107 ^ 

803.4 


138-135Ba 

39(4) 

0.76 

0.072 

-0.000161 ^ 

96.1 


138-134Ba 

114.1(1.0) 

0.76 

0.040 

-0.000217 ^ 

-386.7 


136-134Ba 

-84.8(8) 

0.76 

0.005 

-0.000110 ^ 

805.6 


137-135Ba 

-75.3(5) 

0.76 

0.017 

-0.000108 ^ 

813.8 

6s5d ^D 2 

— 6s6p 







138-137Ba 

69(3) 

0.74 

0.055 

-0.000053 ^ 

-536.9 


138-136Ba 

- 63(1) 

0.74 

0.035 

-0.000107 ^ 

831.6 


138-135Ba 

/ 






138-134Ba 

-143(1) 

0.74 

0.040 

-0.000217 ^ 

795.8 


136-134Ba 

- 80(1) 

0.74 

0.005 

-0.000110 ^ 

761.1 


In theory, the AA^ms parameter is isotope-independent and its value should be identical 
for a given transition. Table IIXI clearly reveals the incompatibility of the AiTMS values 
within a given transition. 


CONCLUSION 


The analysis of the 6s^ ^S'q— 6s6p 3,ipo 0^2 3p^ transitions, gave us conhdence 

in our calculations of the electronic factors. At the end of their paper, Dammalapati et al. 191 
suggest that “the nuclear spin gives rise to an additional contribution to the IS” for odd 
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isotopes, but on the basis of our new results — and also the presence of discrepancies for 
even isotopes — their statement is open to doubt. However, complementary investigations 
would be highly valuable. For instance, with respect to the large MS effects found for these 
transitions, further investigation of the convergence of the ab initio parameters would be 
welcome to conhrm our predictions. The convergence of the electronic parameters with 
respect to more elaborate correlation models should be investigated to rehne the estimation 
of the accuracy of our electronic parameters. Unfortunately, the present calculations have 
reached the limits of the current computational limits. Furthermore, by studying other 
transitions, more conhdence could be obtained on the nuclear mean-square radii between two 
isotopes. To continue this study of the nuclear mean-squares charge radii of barium, many 
investigations remain possible. One of the other possibility would be to reinvestigate the 


46, 


47l | and only few relativistic 


Ba II system, that also presents many experimental studies 
calculations 21]. Another interesting track would be to (re)investigate experimentally isotope 
shifts of transitions involving the 6s5d ^Dj levels. 
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